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Abstract 

P i' 

IjT^' We further investigate, in the planar limit of A'^ = 4 supersymmetric Yang Mills theories, the 

(-H ' high energy Regge behavior of six-point MHV scattering amplitudes. In particular, for the new 

Regge cut contribution found in our previous paper, we compute in the leading logarithmic ap- 
proximation (LLA) the energy spectrum of the BFKL equation in the color octet channel, and we 
calculate explicitly the two loop corrections to the discontinuities of the amplitudes for the tran- 
sitions 2-^4 and 3^3. We find an explicit solution of the BFKL equation for the octet channel 
for arbitrary momentum transfers and investigate the intercepts of the Regge singularities in this 
QQ , channel. As an important result we find that the universal coUinear and infrared singularities of 

f^ ■ the EDS formula are not affected by this Regge-cut contribution. Any improvement of the EDS 

formula should reproduce this cut to all orders in the coupling. 



(N 
> 

0^ 



O. 

^ ; 1 Introduction 

In a recent work [1] we have investigated the high energy Regge behavior of MHV scattering amphtudes 

k> , in the planar hmit of iV = 4 supersymmetric Yang MiUs Theories, and wc have found that, for n-point 

}_( ' amplitudes with n > 5 beyond the one loop approximation, the simple factorizing structure of the 

5t ! Bern-Dixon-Smirnov (BDS) conjecture [2] is not valid. In detail, it was shown that for the cases of the 

transitions 2-^4 and 3^3 their factorized form is violated by Regge cut contributions which satisfy 

the BFKL equation [3] in the color octet channel. These terms are obtained from specific single energy 

discontinuities, and in the scattering amplitudes they become visible in particular physical kinematic 

regions only. In the one loop approximation, these terms are correctly contained in the BDS formula, 

but in higher orders they cannot be cast into the simple exponential form conjectured by Bern et al. 

In this paper we further investigate these Regge cut contributions. We study the BFKL equation 
in the color octet state, and we compute the two-loop expressions for the 2^4 and 3 — > 3 amplitudes. 
In particular, we show that the coUinear and infrared divergences of the BDS formula are not affected 
by the Regge cut contributions. 

The paper is organized as follows. In section 2 we briefly review the derivation of the factorization- 
breaking contributions, and we write down the expression for the Regge-cut contribution, using the 
calculus of complex momenta. Sections 3-5 are devoted to the detailed investigation of this Regge-cut 
contribution: we first (section 3) study the structure of the infrared singularities, we then (section 
4) compute the two loop expressions for the cut contributions, and finally in section 5 we obtain 
the explicit solution of the BFKL equation for the octet channel. In the final section we present 



conclusions and further strategies. Solutions of the BFKL equation for the forward case are presented 
in an appendix. 

2 The Regge cut contribution: review and representation in 
terms of complex momenta 

In our previous paper we have studied the high energy Regge behavior of scattering amplitudes of A^ = 
4 supersynimetric Yang Mills theories. In the leading logarithmic approximation (LLA) we can make 
use of the QCD calculations since the supersynimetric partners of quarks and gluons do not contribute 
(in this limit t channel exchanges with the highest spin dominate) . We now summarize the main results 
of [1] . For n-point amplitudes with n > A, the high energy scattering amplitudes can be written as sums 
of separate pieces (named 'analytic representation' or 'dispersion representation'). This decomposition 
reflects the analytic structure required by the Steinmann relations [4] . The different terms appearing in 
this representation can be computed from single energy discontinuities ('imaginary parts') or multiple 
energy discontinuities. To be definite, we consider the 2 —^ A and the 3 — > 3 scattering amplitudes, 
illustrated in Figs. [1] and [21 For the 2^4 scattering we are interested in the kinematic limit (double 




Figure 1: Kinematics of the 2^4 amplitude 
Regge limit) 

S > Sl,S2,S3 > il,t2,t3, (1) 

whereas the 3^3 scattering process will be studied in the limit 

S > Si3,S02 > Sl,S3,t2 = ipA - PA' - ^2)^ > ^1,^2,^3- (2) 

The analytic decompositions are illustrated in Figs. [3] and HI 

In the physical region where all energies are positive, there are, both for the 2^4 and for the 
3^3 case, substantial cancellations of the Regge cut contributions between these five terms: in the 
sum their imaginary parts cancel, and the amplitude takes the well-known factorized Regge form. 
In other physical regions, however, where some energies are positive and others are negative, the 
cancellations are less complete, and pieces become visible which do not show up in the region of only 
positive energies. For the 2 — + 4 case, the physical region of interest is 

S > 0, S2 > 0, Si < 0, S3 < 0, Soi2 < , Si23 < 0. (3) 

Here non- vanishing discontinuities are only in s and S2 , and both of them contain a new term which 
violates the simple factorizing form. It contains a Regge-cut structure which is described in terms of 
the color octet BFKL equation. We illustrate these discontinuities in Fig. [51 




Figure 2: Kinematics of the 3 — > 3 amplitude 




Figure 3: Analytic representation of the amplitude Ah^i 
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Figure 4: Analytic representation of the amplitude M^^^ 
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Figure 5: (a) the S2 discontinuity for the 2 
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Figure 6: The S2 discontinuity for the 2 — > 4 amplitude. The big blob denotes the BFKL gluon Green's 
function in colour octet state including the reggeization of the lines with momentum fc, g2 — k, k' and 
q2 — k'. The dots indicate the production vertices in eqs. (O and P^ . This is equivalent to Fig. [5] (a) 
indicating the momentum notation. 



For the 3 — > 3 scattering amplitude the corresponding region is 



s > 0, 4 > 0, si < 0, S3 < 0, si3 < , so2 < 0. 



(4) 



The non- vanishing discontinuities belong to s and <2 > 0, and they, again, contain the Regge cut 
pieces. 

In [1] we have compared these results with the expression given by Bern et al. Whereas for the 
2^2 and 2^-3 amplitudes the QCD results are in full agreement with the EDS formula, the 2^4 
and 3 — > 3 EDS amplitudes are correct only in the one loop approximation. For two or more loops, 
the Regge cut piece cannot be reproduced by the EDS expression. As explained above, this implies 
that the EDS formula (in LLA) still gives the correct result in the physical region where all energies 
are positive, but it fails (beyond one loop) in the regions ^ and (|4]). 

In the following we shall investigate these Regge cut pieces in more detail. Rather than returning to 
the five terms illustrated in Figs. 3 and 4, we directly present an explicit Feynman diagram calculation 
of the single energy discontinuities in S2 and s for the 2 — > 4 amplitude and in t'2 and s for the 3 — > 3 
amplitude. Let us begin with the S2 discontinuity illustrated in Fig. [Sj Here the blob in the center 
denotes the EFKL Green's function in the color octet channel which sums the s-channel emissions in 
the center of Fig. [5^, and on both sides we have to convolute this Green's function with the 'impact 
factors' <&! and $2. Introducing complex momenta 



r^X ~f ^"J-i/i "^ — '^X lif\iij 



(5) 



and making use of the expression for the vertex describing the production of a gluon with definite 
helicity (cf. eq.(6) of [1]): 



C^(g2,gi)e*(fci)^x/2M1 
we obtain for the production vertex to the left of the Green's function 



%/2. 



gi(g2 ~k)* 
(k + ki)* ■ 



(6) 



(7) 



Here we have used that, in Fig. [SI the gluon with momentum k + ki is on shell (we consider the 
discontinuity in S2), and at the upper vertex where the gluon with momentum k is attached the 
outgoing gluon helicity is conserved. Since the scattering amplitude T2^„ for the case of the maximal 
helicity violation (MHV) can be written as [2] 



T, 
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Born . ^. 



(8) 



k+kn 



q.^'VW'^i— q2-k q^WvUci;^ 



q^-k 




Figure 7: Decomposition of the production vertex in eq.Q. The dots denote the effective production 
vertex in eq.®. 



we will, throughout our paper, consider the factor M2^„ only. We, therefore, separate the production 
vertex of the Born approximation and rewrite ([7]) as: 



^/2^<i>l(fc,q„9l) 



(9) 



with the impact factor: 



^i{k,q2,Qi) 



kl{q2-ky 



= 1 



k*ql 



qlik + kiY 



(10) 



In the following we shall work with this impact factor. 

In order to make contact with [1] we should note that, with the result in the second line of eq. H 
the production vertex in ([7]) can be written as a sum of two terms of the form 



^gi(q^-fcr ^ ^'Zi'Z2 



(fc + fci)* 



kl (fci + k) 



ql ^{k, + k)k* 



K 



(11) 



We illustrate this structure in Fig. [T] The first term is 'local', i.e. it has no further dependence on 
the internal momenta, whereas the second one is 'nonlocal'. 

Similarly, on the right side of the Green's function in Fig. [6] we have 



V2. 



fc2 — k' k2 



(12) 



with the impact factor 



*2(fc', 92^93) 



= 1- 



k2{k'~q2) 
q2(k'-k2) 
k'qs 



{k' - k2)q2 



and 



V2 



(g2 - k')ql 
ki-k' 



V2 



<i2% 



ilk' _^(fc2-fc')* 



(13) 



(14) 



V2 - «;■ k2 (fc2 - A:')2 ' k2 

The discontinuity in S2, to an arbitrary loop accuracy, of the amplitude M2_>4 in LLA then has 
the form: 

TT ''' ^^^ Ja-ioo 27ri V/i2 

where the i2-channel partial wave is 



f2it0), 



(15) 



/2(a;) ^a.qlj d'-''kd^~''k' $i(fc, q,,q,) Gl'-\k, k' , q,) $2(fc', 92, ^s) 



(16) 



and 



Oif = ,^ sn o, ,a=— — (47re 'j 



(27r)2-2e '"' 27r ^^"" ^ ' ^^''^ 

The overall factor q| in front of the integral in p6p takes into account that the Born approximation 
of the amplitude contains the pole l/|g2p but in our calculations we are interested in the scattering 
amplitude M2^4 with the Born factor being removed (cf. [8|). The Green's function Gli {h^h! ^q^ 
satisfies the BFKL equation for the color octet channel (putting e = 0): 



k (fc-qa) fc (fc-92) ^ ^ 



(18) 



where K'^^^^ denotes the BFKL kernel in the color octet channel, containing both real emission and 
the giuon trajectory, and the convolution symbol stands for ^ = J jf-w- Using complex momenta 
the kernel can be written in the form: 

K^'^\k,k';q,) = S^^Hk k') (c.(-|fcP) +^(-|g2 - kf)) + f^'' - ' |fe y ^* ^"" ' (19) 



where the gluon trajectory is 



, 1 k"^ 

^^--'i-ln^ 



Lo{—k ) = a 



(20) 



In contrast to the color singlet BFKL kernel, the color octet kernel is not infrared finite and needs 
to be dimensionally regularized. It is convenient to separate the singular pieces by writing the octet 
kernel as 



K'^^^Hk,k';q.2) = 6(^\k~k') 



'hql) + 1 H-k') + ^i-{Q2 - kf) - M-qD) 



lK<^'\k,k';q,) 



2 

6^^Hk-k' 



' I 2n «, fe^(92-fe)^"' 



;K('\k,k';q,) 



(21) 



In this expression, K^^'>{k, k'; Qj) denotes the color singlet BFKL kernel, and infrared singularities are 
contained in the trajectory function w(— q^)- 

Inserting this form of the octet kernel into (|15p , the discontinuity takes the form: 



27ri \ij,' 



-%,M,^,^sf'^^ I ;E?(^) /2M 



(22) 



where the reduced partial wave f2{^) is given by 

h{u;)^a,ql [ d'-'^k d'-^'k' 'P,{k,q^,q,) G^{k,k' ,q^) ^^{k' ^q^^qs) 



(23) 



The Green's function G^{k,k\q2) satisfies the BFKL equation ([18]) with the reduced kernel 
k^k.k'.q^) = -5^^\k-k')Un^^^^ + \K<^^\k,k';q,). 

^ 9292 ^ 



(24) 



From the explicit form of the function /2(w) and of the impact factors $i one sees that there are 
potential divergences only for |A;| ^ |fe'| -^ (and not for 1^2 — k\ ^ \q2 — k'\ -^ 0). The one loop 



contribution to the partial wave, /j , takes the form 

In our previous paper [1] we isolated the term which violates the BDS factorization ansatz. This 
term, named Vcuti is contained in (J16p : in the impact factors $1 and $2 one retains only the nonlocal 
pieces (cf. (|lip and Fig. [7]), and one subtracts the Regge pole contribution. The one loop contribution 
was given in eqs. (94) and (95) of [1]. In the normalization of (|16p it reads: 

k*ql 1 kqs ^ a / ^ gjgj 

q^ik + ki)* k^{q^-kY q2{k-k2) 2 ^ "" (fci + fea) V^ 



^^<li d^^'-k—^^i^— -^——- ...,.. - :^ In „. ^T..... --] (26) 



and it was shown to coincide (apart from an overall factor) with the phase factor C in eq.(75) of [1]. 
In this paper we address the discontinuity in $2 , for which we do not need to split the impact factors 
into local and nonlocal pieces. However, for the discontinuity in s we will come back to the result 

dMD. 

In the following sections we will study the reduced partial wave (|23|) and the reduced kernel ([24|) 
in some detail. First we will investigate the infrared properties and show that the infrared divergence 
is contained only in the one- loop approximation (PS)) . i.e. the reduced kernel is infrared finite and 
introduces no further divergences. This implies that the divergent term ~ 1/e is not renormalized. 
We will then compute explicitly the two loop approximation to the reduced partial wave. Finally, we 
will return to the reduced kernel and calculate its eigenfunctions and eigenvalues. 

3 Infrared properties and eigenvalues of the octet kernel 

In this section we concentrate on the infrared properties of the octet kernel. By investigating the 
most singular part of the reduced partial wave ([^5]) . we find the exact expression for the eigenvalues, 
and we prove that the infrared singularity of the reduced partial wave coincides with the 1/e pole of 
the one loop approximation. The exact solution of the color octet BFKL equation will be derived in 
section 5, and it allows to find a closed expression for the reduced partial wave. 

The starting point of our further discussion is eg. ((23)) . During this section we will denote the 
transverse momenta by p,p' and p — Px + iPy,P* — Px — ipy The Green's function G^j satisfies the 
'renormalized' equation 

^G^{p,p',q) = ^- -^ S^{p - p') - aH G^{p,p' ,q2) , (27) 

p [q-p) 



g = ln '^",i|;^' + J"^rr P(g-Pr + ^^^P*(g-p) + 27. (28) 



where 

bp |g-pp , 1/p ln|pp , l/pMn|pp 

\q\^ q*-p* 2 q-p 2 

As we stated before, the most interesting region is the infrared divergent region \p\ ~ \p'\ <C \q\. In 
this asymmetric kinematic it is possible to find eigenvalues and eigenfunctions of the reduced kernel. 
First, the expression ((23|) for /a is simplified: 



J2{uj) ^a, j d'pd'p'gUpJ), (29) 

and 5ij satisfies the equation 

^9uj{p,l^) = ■p-f^(5^(p-p') -aHg^{p,p) . (30) 

\P\ 



Here 

H = ]n l«|2 + 

p 2 ^ p* 2 



H^ln\p\^+'-'^P+V^P* + 2, (31) 



and 7 — — V'(l) is the Euler constant. 

The Hamiltonian for the octet quantum numbers has the property of the holomorphic separabihty 



Inp 1 Inp 



H^hs + K, hs = lnp+^ + -^p + j. (32) 



The holomorphic Hamiltonian /ig is slightly different from the corresponding Hamiltonian for the 
singlet case 

h = -— = Inp + - (In p) p + 7 . (33) 

2 p 

The difference is also in the normalization conditions for the wave functions in these two cases 

m\i^ Jd'p^*\p\'^, wnipKL^ Jd^p^*\p\'^- (34) 

The eigenfunctions belonging to the principal series of the unitary representations of the Mobius group 
in the holomorphic subspace have the different form 

The eigenvalue of the total Hamiltonian for the octet case is given by 

To verify this result we act on the wave function f{k) with amputated propagators with the Hamil- 
tonian regularized by a mass parameter p,^ 

1^1 \u\2iv ( ^ \ f dk kk* + k*k ,,/,2wfk \'^ 



^^ ^"^ fi^ '^' " Uv J 27r|/fc'|2 (|fc-fc'|2 + ^2) l^'l " [k'* 

= In L^ |A:|2- f A^ ^ ^ f\x I' '^^^' ^^ " ''^""^^ {l-^^+ f) (^A:'* + fc*fc') fc'" 



Ai^ ' ' V^V 7o 7 27r(|fc'-x/c|2+x(l-a;)|fc2|+a;;u2)2-»^+f 



//2 ,,2 I n^ / ■? \ l-'ii'+'f ' ' \ J.* 



4 •"' [^-^ + w 

One immediately sees that the result is finite when ^^ is taken to zero. 
From this expression we can obtain the eigenvalue 

Eun = U^-^ ^-)+^(l + it.+ |n|/2) + ^(l-zt.+ |n|/2)-2^(l) 

= ^xP{l + iiy + n/2) + 3iip{l + iiy-n/2)-2tlj{l). (38) 

In particular, for n = 0, 1, we have 

E,,o = 2ntP{l + iiy)-2tl^{l), 

E.i = i^^ + 25ft^(i + zj.)- 27^(1). (39) 

2 j/^ + -. 2 



LU 




Figure 8: Spectrum of the lowest eigenvalues 

The corresponding lowest energies are i?'-*'-' = and E^^^^ =2 — 41n2<0. In Fig. [5] we show the in- 
dependent eigenvalues for different values of n. Thus, the ground state energy corresponds to \n\ = 1, 
as it was in the case of the colorless Odderon state [5] . In section 5 we will re-derive this spectrum by 
solving the eigenvalue problem exactly. 

The solution of the equation for the Green's function can be found with the use of the completeness 
condition for the eigenfunctions 



9uj{p-,'0) 



1 



1 



2U/|2 



\P?\P 



.(0-0') / ^ 



dv 

2^ 



HI 



1 



- ijj{v, n) 



(40) 



Here 



ijj{v,n) = — - 



t^vn — 



■^■K^ 



87r2 



(2V'(l)-3fi?A 



- ni' 1 + iin - 



(41) 



Inserting this expression into (I29p one sees that /2, in this approximation, reduces to the 1/e-pole of 
the one-loop expression: the reduced color octet BFKL Green's function only leads to finite correc- 
tions to the one-loop result and introduces no further infrared divergences. Therefore, the divergent 
contribution is 

1- -^ i ^ ij(t2)+w„(o,o) _ _a 

e ^ ^ e ■" 



-%,M2^. 



IT 



2^i\div 



J{t2) 



(42) 



where u}{0, 0) is the leading eigenvalue of the 'reduced' octet BFKL kernel discussed before (note that, 
in this 'infrared' approximation, the impact factors are equal to 1, and the solution belonging to the 
eigenvalue a;(0, 1) does not contribute). Thus, the divergent contribution ~ 1/e is not renormalized: 
its asymptotic behavior corresponds to the usual gluon Regge pole. The reason for this is that the 
collinear and infrared divergences are factorized and that the EDS ansatz is valid in the one-loop 
approximation. 



4 Contributions at two loops in the perturbative expansion 

Let us return to the reduced partial wave f 2(1-0) in eq. (|23p and compute the first terms in the per- 
turbative expansion. The one loop approximation has already been given in (|25p . As mentioned 
before, it contains the infrared singularity coming from the region |fc| ^ 0. In this section we iterate 



10 



the integral equation for the Green's function (|27p inside the partial wave and, using the calculus of 
complex momenta, compute the two loop expression. Starting, in Fig. [6l from the impact factor on 
the right, given in eq. ()13|) . we have the following expression for the first iteration (using the reduced 
color octet BFKL Hamiltonian in eg. ([24]) ): 



192 -fc P 
192 Pa^^ 



H$2 = -a In ' W^' ;' $2(fc,g2.g3) 



[ d^k' kk'*{q* - fc*)(g2 ^ k') + rfc'(g2 - kM - k'*) 

^ ""J ^^ (|fc-fc'P+^2)|fc,|2|5^_fc,|2 'J'^lfe , 92, 93) ■ (43) 

Here /z^ plays the role of an intermediate infrared cut-off which will be removed at the end of our 
calculations. The result of integration can be written in the form 

H^2 = -lx{k), (44) 

92 ^ 

where 

XC^) = -« I ^ In i„ I4IL r^ + TT In ^7^ + 7-71 J:T In ^^ U2 ■ (45) 



Next wc perform the integration over k with the impact factor on the right hand side in eq. ljlOp . and 
we obtain for the two- loop approximation of the imaginary part in the S2-channel in (|22p : 

As. = r^'°°— f-V/2(^) 



n ^ ko , f d k s /, ^ , s 

= -o92 — a lns2 / — ^^ — $1 fe,q2,9i X fe • 46 

2 92 J 2'K\k\^\q2-kY 

With the use of complex number algebra this two-loop expression can be reduced to the form 



. TJ" 2 1 f d^k 

As2=-^a lns2 / -^PW^ 



(47) 



where 



1 


1 ^ 


In'^^" 

I92PI 

1 " 


k\'\k2\' 
k~k2\'' 
\ 1 I92 - 

Mni— 12- 
) In 1'^^ - 


-k2\^\k\^ 
\k^k2\'' 
-k\^\k\^\q2 
19211^^- 


-k2\ 




\k\' 
1 


k{k* + kl)J 
1 1 




fc- 92 
1 


k* k-q2 

1 


1 1 


'\k2? 


k — k2 


k* + kl k 


-k2k* 





(48) 

One can easily verify that the ultraviolet divergences cancel. Also, in agreement with the previous 
sections, the divergence at fc = is absent. The above integrals over k can be expressed (with the 
shift A; — > fc -|- c) in terms of the following expression: 

^("'^)"/ .(fc-aKt-6V '^'^'^- ^''^ 

To regularize the ultraviolet divergence we introduce the cut-off 

|fcp<A2, (50) 
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which at the end cancels in the expression for As^ ■ One can then write / in the form 

In^A^ 
/(a,b) = ^- + /,(a,6) (51) 

and use further the regularized value fr because In A is canceled in the final result. To calculate this 
function we take derivatives in the complex coordinates a* and b 

^^--^l^H^|/--^l^|5^ (52) 

After integrating these expressions we obtain 

. , , , f dx f dy 
fr{a,b) = - / b^lna;- / ^Iny 

Jo X- — Jo y- -b 

+ InlalMn ^* +ln|bPln— -l-a^(ab*). (53) 

a* — b* — a 2 

The last term was obtained as an integration constant: it can be determined from the conditions that 
it must depend upon a and fe*, and the full function / should depend on the invariants 

\a\^=aa\ \b\^^bb*, ab* ^ ab - i[a,b]3 , a* 6 = ^-|^ . (54) 

Moreover, from dimensional considerations it follows that it contains the term i In s, where the 
invariant s has the dimension of A^ . 

The function fr can be expressed in terms of the Spence's function (dilogarithm) 

fria,b) = -iz2(^^^-Lz2Q^-lln2(a6*) 

+ ln|a|2ln-^ + ln|6|2ln-^, (55) 

a* — b* b — a 

where 

Li2iz) = - H-^^^^^-^dt. (56) 

Jo ^ 

Note that the above expression for / has the following properties 

/(-a,-b) = /(a,b), r{a,b) = f{b,a). (57) 

In some particular cases it can be simplified. For example, 

We shall use also the values of the integrals 

7r|fc|2 ^ |c|2 "2 "" |c|2 ' y 7r(fc-a)(fc-6*) ^ ""la-bp- ^^^' 

With these results we can calculate the two-loop contribution to the imaginary part of the ampli- 
tude in the S2-channel: 



'' - lnf^lnLlL-/,(-g2,-q2-fei) + /,(-fc2,-fci-fc2 



a2 1ns2 TT |(j2p |g2p 



"1^ r-[2 1^ r-[2 +/r(Q2,0-./r 92,-^1 +/r(g2 -^2,-^1 - fc2 

92 92 r 



k2? ^^ |92-fc2p |fc2 

ki+k2\^ ^ \q2\ 

fr{k2 - 92, -fel - 92) - /r(fc2 - 92, -(I2) - 2/r(0, -fcl - ^2) + 2/,(0, -^2) . (60) 



-fr{q2 - fc2,-fc2)+ln ,,' ,'', ,2 In '"" ,„,',' " +./.(fc2,-fci) - /.(fe2,0) 
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Using the following properties of dilogarithms 

'1\ .. . . 1 



i?;2 (1 - z) = - Li2{z) - ln{l - z) \n z + (^2 , 

LiafY^j = -Li2iz) - -\n^{l- z), 

Li2 (y3^) = ^^2(z) + ln(l - z) H-z) - i ln2(l - z) + (2 , 

ii2 ( ^Y^ ) = Li2{z)+Hl-z)lnz-^ln^z-(:2, (61) 



we can simplify the following sums entering in 2As 



- friq2,Q2 + ki)- friq2,-ki) = Inj/jip ln|9ip-2C2, 
fr{k2,ki+k2) + frik2,-ki) = - hi |fci ^ In |gi - (73 P + 2^ 
fAq3,-ki-k2) + friq3,qi) = - In \qi\^ In \qi - qsl'^ + 2C. 



2 ■ 



2 



-friq3,-k2)-friq3,q2) = ln|A;2|Mn|92p-2C2. (62) 

The final result for Ag^ can be written in the very simple form 

. a% 1 Igl-g3p|g2p , |gl -g3p|g2p .„„. 

As^ = -TT— In S2 In I |2|, 12 In , .ai, 12 • (63 

It is symmetric with respect to the simultaneous substitutions 

fci ^ ^2 , Qi ^^ —93 • (64) 

In a similar way we can calculate the discontinuity in the s-channel. Starting, in Fig. [5)3, from the 
gluon ladders in the ii and the ^3 channels, we invoke the bootstrap equation. This equation allows 
us to write, instead of the gluon ladders, simple Regge pole exchanges. The resulting s-discontinuity 
has the same form as the S2 discontinuity with the impact factors <&! and $2 being replaced by 

^' k* + kl q* ' ^' 92 k'-k2- ^ ' 

One easily verifies that these modified impact factors coincide with the nonlocal pieces of $1 and $2 
in eqs. pH)) and P^ . We also note that Ag can be obtained from As^ by substituting 

k'l ^ -gi , ^2 ^ g3, (66) 

and by changing, inside Fig. 6, the integration variables k —> q2 — k^ k' ^> q2 — k' . In fact, in the two 
loop approximation, Ag coincides with Ag^ '■ 

1 Igl^g3p|g2p i(,^^ 

g3 r Fl 



(68) 









Ag^ 
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1 , lgi ~ 
in S2 In — ; — 

|gi 


12 

gal 

P|fc2 


|g2p 


due to the 
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-momentum conservation 
















fcl + fc2 = 
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Figure 9: The t'2 discontinuity for the 3^3 amphtude. 



Let us now consider the non-factorisable contribution for the scattering amphtude 3^3 (Fig. [2]). 
In this case we, again, have the imaginary parts in t'2 and s-channel. For the imaginary part in the t'2- 
channel (Fig. [9]) we have, on the left side, a shghtly modified impact factor, $1: the two corresponding 
impact factors are 

k* kl ^ ^,^ k' -q2 fc2 



*i(fc, 92.91) = -- 



91 92 



^, ^2{k) 



k' — fc2 52 



(69) 



and, therefore, the infrared divergence at A: = is absent. For completeness, we first list the one loop 
results. The one loop result for the partial wave (analogous to (|^5|) ) is: 



^1- 



a^qi / d^ ^'k $i(fc,q2,gi 



1 ^ n ^ "1 ^2(91 +93 -92)^ 



(70) 



Note that, in contrast to the 2^4 amplitude, there is no infrared divergence. Next we quote the term 
analogous to ((26|) which is obtained by retaining, in the impact factors, only the 'nonlocal' pieces: 



- 2 / j2-2e, 9i(92 -k) 
a^q2 / d k 



1 



fc93 ^ a , 92(91 +93-92)^ 

(72*(9i-^)*^'(92-fc)2g2(fc-92) 2 "" klkl 



(71) 



It coincides (up to an overall factor) with the exponent of C" in [1] (eq.(80)). In contrast to the 2^4 
amplitude, this one loop result, again, is infrared finite. 

For the two loop calculation we proceed in the same way as for the 2 — > 4 case. Using our previous 
results for the function x(fc) we obtain 



Ah ^qja \nt'2 



d^k 



— $i(fc,q'2,gi)x(fc). 



2 ' 27r|fcP|92-A:|2 92 
With the use of complex number algebra it is possible to transform this expression into the form 



(72) 



A,, 



-a^ Inii 



d2fc 

27r 



?(fe), 



(73) 



where 



p(fc) = 

+ 



|g2-fcp|fc2p 

k{k*-ql) k{k*-qX))^^^\q2\^\k-k2\^ 



In 



1 



1 



1 



2|l,|2 



In 



l92-fc2p|fe| 

|fc-92P fc-92fc*-9ry'"l92P|A:-A:2|2 

^ 1 ' ' -m. 



k-k2k*-q*2 k-k2k*-q*J |fc|2|fc-g2p|A:2|2|g2-fc2p 



92|1fc-fc2|^ 



(74) 
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The integral over k can be expressed in terms of the function /r(a, b) introduced above: 



2 . . . , , , ^ . . , ^ , 192^ , \k: 
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A' = /r(-fe2,-A;2 + qi) -/r(-q2.-g2 + 9i) -In-r-r? 1^ I — I? +/'-(-g2.0) 



9 2 

-/r(-fc2, gs) + /r-(92 - ^2 , -^2 + Ql) - /r(92> +9l) + 1^ 4 ^^ tI 

92 '^l 

MnM|iin(^l^_2/.(0,,,-fc2) 

92 93 

+ 2/r(0, 92 - ^2) + /r(fc2, 9l) - /r(fc2, 92) + /r(fc2 - 92 > -^2 + Ql) - /r(fe2 - 92: 0) • (^5) 

With the use of the identities for the sums of the functions /(a, b) hstcd in ([5^ . we can significantly 
simplify A^i : 

. 7^ 2 1 ,/ , |g2-gi-g3p|92p , [92 -91 ~93P|92P .„„. 

^*2=~7^ lnt2 In In ■ — — — . 76 

' 4 |fcir|K2r l93rl9ir 

Thus, Afi is different from A^^ and Ag by the substitution Qi <-^ — fei. In fact, one can also verify 
that the same result is obtained for the imaginary part in s for the 3^3 transitions. 

^^3 ^ „- 2 1^,, 1^ ''^^'f.^f'"''^^'' In ''^^','^\~'^^'''^^'' . (77) 

4 |A:ip|A:2p l93p|9ip 

As indicated before, all the two loop results are infrared finite and, hence, do not affect the infrared 
singularities in the BDS formula. In the next section we find the explicit solution at all loops. 

5 Solution of the BFKL equation in the octet channel 

In this section we solve the eigenvalue problem for the reduced color octet kernel and derive all-order 
expression for the 2 — > 4 and 3-^3 amplitudes in the leading-log approximation. For the eigenvalue 
problem it is convenient to return to the symmetric notations of the momenta pi — p, p2 = q — p and 
write the homogeneous BFKL equation for the wave function / with the removed propagators in the 
octet channel as follows 

Ef{pi,P2)=H,f{pi,p2), (78) 

where H has the holomorphic separability property 

H = h + h* , h — In ^ 7: { Pi In P12 1- P2 In 012 — ) + 7 • (79) 

92V Pi P2J 

With the use of the relations (see [6]) 

\n{z^d) = In z + i {^/j{zd) + ^{-zd + 1)) , \n{d) = - Inz + ^ (V^zO + 1) + ip{-zd)) (80) 
one can transform the holomorphic Hamiltonian to the form 

h=^\nq+-{\n {pipi2) + In (P2/O12)) + 7 • (81) 

By introducing the conjugated variables 

Pi « ^ -^2 /oo^ 

y = — , 9 = — = -^ — P12 , (82 

P2 oy q 



h can be simplified as follows 



h^^ (Hy^d) +lnd)+j^^ {^{yd) + ^(y^ + 1)) + 7, (83) 
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where we neglected pure imaginary terms which cancel in H . 

Thus, the solution of the homogeneous BFKL equation in the momentum space can be found in 
the form 

Un{k,q) = ' 



q — k 
The corresponding energies were calculated above 



q* — k* 



E^ 



^ ( «^ + 2 ) + "^^ 



The orthogonality condition for the above wave functions is 

(fik 



7r|fc|2|g~fc|2 

Their completeness condition can be written as follows 



ve functions is 

rt'n'i^^ <?) Iun{k, q) = 2ti6{v' - v) 5n',n 



2V(1) 



E 



dy rUk'. <f) fun{k, q) = 2^'5'{k' - k) 



\k?\\l-k\ 



(84) 



(85) 



(86) 



(87) 



Therefore the Green's function for the i-channel partial waves is 



1 



Gu,{k,k'\q) = —ir 



where 



27r2 |fc|2||g-fc|2 



uj^v, n) 



E 



n— — oo 



g^Na 



dv 



f:n{k'.<f)Un{k,q) 

uj — uj{i', n) 



8^2 ^-■ 



(89) 



With these results we can find explicit expressions for the S2-discontinuity of the 2^4 scattering 
amplitude and for the ij-discontinuity of the 3^3 scattering amplitude. Starting from eq. ipS)) . we 
have to convolute the octet channel Green's function with the corresponding impact factors. Returning 
to Fig. [5] and to the notation of section 2 we have to calculate the integral 



X2 



1921 



<fk' 



2^ |fc'|2|g2 



q-i-k' 



l\ iv+^_ 



92 



L/*\ "■"- 



k2{k'-q-2) 
(k'-k2)q2' 



(90) 



The simplest way to calculate X2 is its differentiation in /cj with the subsequent integration, which 
gives 



^^ 2(^.-|)Vfc2 

In a similar way the integral over k gives 



1 1 /93\ lis 



Xi 



2(i;/+f) V fci 



91 



kl 



-iv+^ 



(91) 



(92) 



As a result, the imaginary part of the production amplitude in S2 for the transition 2^4 takes the 
form 



-^s.M2^, = ^4'^^ £ {-IT Regs 



dv 



qlki 



" \ h*n* 






L!{iy,n) 



(93) 



where the regularization refers to the divergence at i/ = 0, n = 0. which appears only in in the one loop 
approximation. In appendix B we compute the one and two loop results (obtained from expanding 



j{i'^n) 



= 1 + hiS2Lo{iy,n)), and verify the agreement with ([^5]) and (|63p : 



-5., M2-.4 = ? s^^*^^ (In- 



\ki\^\k2 



k2\^^l 



2, ,2 



In S2 In 



|A;i+fc2|2|<Z2|2 \k,+k2\^\q2 



|fc2Pkl| 



In- 



IfciPkal 



(94) 
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In an analogous way we compute the discontinuity in s. In (|90[) we replace the impact factor $2 by 
<I>2 (and similarly for $1 in (|92p ). and proceed in the same way as before. The result can be written 
in the form 



1 °° PC 



dv 



Q3K 

^2 + ^ V^2'?l 



7 \ jl/'+T 



j{v,n) 



(95) 



with the regularization prescription Regs for the singularity at ly — 0, n = which, again, applies to 
the one loop approximation and takes care of the difference between the discontinuities in S2 and s. 

As a result, the production amplitude 2 — > 4 in the multi-Regge kinematics with s, S2 > and 
Si, S3 < in the leading approximation can be written as follows 



A2^4 



aBDS 



(l + iA2^4), 



where 



i2-^4 



? E (-1)" 



di> 



MM 

^2 + ^ \k*ql 



^291 



(s^^"^'") - 1) , 



(96) 



(97) 



has no infrared singularities. We mention that in the region s, S2 < and si,S3 > the scattering 
amplitude has the similar form 



A 



2^4 



A 



BBS 



(1-ZA2_4). 



(98) 



We emphasize that the correction A2^4 does not contribute outside these physical regions. 




Figure 10: Dual variables for ^2^4. 

It is noteworthy that if we perform the duality transformation shown in Fig. \W\ (cf. Ref. [7]) 

qi -^ zoi , ki -^ zo'i , qs -^ zo2 , ^2 ^ Z20' , (99) 

and introduce 'coordinate' variables Zi, we see that our results for the imaginary parts depend on the 
anharmonic ratio 

x^^^Hfia. (100) 

Z0'2Z01 

The reason why the BFKL equation in the octet channel can be solved is its invariance under Mobius 
transformations in these z, variables. 
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It is interesting to note that the correction to the BDS formula in our kinematics can be written 
in terms of four dimensional anharmonic ratios [8,9]. In particular, in second order of perturbation 
theory we can write 






-2i7r— In S2 In — -—j- — r^ — In — -—j- — r^ — 

4 \k2r\q1r \kirmr 

a^ , ,, (1-$), (1-$) 

—Li2(\ - $) In -^ '- In -5^ '- + . . . 

4 ^ ' $2 $1 



(101) 



where the dots indicate corrections beyond the leading logarithmic accuracy, and we have used the 
notation 



$ 



SSi 



S012S123 



$1 



51*3 ^ 53*1 

^2 



5012*2 



Sl23*2 



(102) 



An analogous result holds for the 3^3 amplitudes (for details see Appendix B). The discontinuity 
in ^2 of the scattering amplitude 3 — > 3 in the multi-Regge kinematics with s, ^2 > and si, S3 < in 
the leading approximation is given by 



-K 2 ^^ 47r 2 



E (-l)"^e.g. 



dv 



(<m 



,,2 I " \ h*]c* 



1^]" %M-.«)^ (103) 

/C2/C1 



(where, in this case, the regularized integral over v for n = and a ~ Q does not contain any 1/e 
divergence), and the 3 — > 3 amplitude takes the form 



where 



A3^3 = f E i-^y 



A3^3 = ^^^3^(l + *A3^3): 



dv 



( gsgi \ ' Aagi 

,2 + 1^ U2*fclV U2fcl 



{t 



luj{u,rL) 



(104) 
(105) 



In the region s, <2 < f^^id si, S3 > we can write 



A3^3 = ^rir3''(l-*A3^3). 



(106) 




Figure 11: Dual variables for A 



3^3- 
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Similarly to the 2^4 scattering amplitude, if we perform the duality transformation (see Fig. fTTj) 

qi -^ ziQ> , fci -^ zoi , 93 -^ zo2 , k2 -^ Z20' , (107) 

the imaginary parts of the 3 — > 3 scattering amplitude in the <2 ^nd s channels depend on the same 
anharmonic ratio (|100[) . Again, the corrections to the BDS formula can be expressed in terms of four 
dimensional anharmonic ratios. 

From these results for the 2^4 and for the 3^3 amplitudes we conclude that the infrared 
structure of the inelastic amplitudes is given correctly by the BDS expression, whereas the finite 
factors are correct only in the one loop approximation. 

6 Conclusions 

In this paper we studied, in the leading logarithmic approximation, the cut contribution which, in our 
previous paper, was found to violate the simple Regge factorization of the BDS formula. As a main 
result we have verified that the factorization of universal infrared singularities is not affected, i.e. the 
violation of the BDS formula is in the finite pieces. We have computed the energy spectrum of the color 
octet BFKL Hamiltonian, and we have concluded that the infrared divergent gluon trajectory can be 
separated from the finite remainder of the BFKL Green's function. We have explicitly computed the 
two loop approximation of the Regge cut piece. The integral equation for the wave function of two 
reggeized gluons in the octet channel is solved explicitly and the intercepts of the Regge singularities 
are calculated. 
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A Forward scattering with color octet exchange 

In this appendix we relate the octet equation in the forward case to the BFKL equation for the singlet 
case. As it will be shown, its physical solution exists only if the gluon Regge trajectory is not expanded 
in a power series on e. 

The relevant integral equation to be satisfied by the gluon Green's function in the octet case reads 

ujG^{p,p',q) = —rj r^S^~^''{p^p')~aH,G^{p,p',q2), (A.l) 

\p\ \q-p\ 

where 

a=^{Ane~^)\ (A.2) 

and 7 is the Euler's constant. From now on we shall use the notation pi = p and p2 = q — p for the 
transverse momenta of the two Reggeized gluons. 

The BFKL Hamiltonian for the channel with the octet quantum numbers can be written in operator 
form as 

He = In -^ h In \p2 H P1P2 + — 7, P1P2 + 27 , A.3) 

M e P1P2 2 p\p2 2 

where we have neglected terms of 0{e) and introduced for the two Reggeized gluons the complex 
variables pk, pi (pi2 = Pi ~ P2) and their canonically conjugated momenta pk,pl.. The dependence on 
the divergence 1/e can be removed by the following shift in the parameter uj: 

uj^uj + -, (A.4) 

19 



which leads to the appearance of a Sudakov-type infrared divergent factor in the amphtude M2^. 



4, 

M2^4 ^ ZM2^4 , Z = exp ( ^ In ^) . (A.5) 



M^ , ^ 7 M^ . 7, ^ P^r, I - In - 

We can then work with the renormahzed Hamiltonian H removing the divergent term 



H = H, + -. (A.6) 

e 

It is known that the BFKL equation in the color singlet channel is Mobius invariant in coordinate 
space. Its solutions are 

\PWP20J \PwP20J 

where m, fa are conformal weights 

1 71 1 71 

m=-+w+-,m=-+ii/-- (A.8) 

and pb is the Pomeron coordinate. The expression (|A.7P corresponds to the three-point Green's 
function with non-amputated legs. The scalar product of two of these functions is defined by 

< Em,fh\\Em' ,fh' >= / (fpi(fp2Em^ffi{piQ,p2a)AiA2E*^,~{PW',P20'), (A. 9) 

where Ak are the corresponding Laplace operators. 

A.l The solution for the octet case at ^2 = 

The solution in momentum space for q2 — in the color singlet case can be obtained by using the 
Fourier transform of the singlet solution (|A.7|) , i. e. 

UMp) = / d^Pi2 d?Poe'^''^ E^nMPio, P20) - P""' b*)""' ■ (A.IO) 



It is convenient to introduce the new function (j)u_n{p) as follows 

frn,fa{P)^W^^AP)^ (A-H) 

with normalization 

^.Ap) = (^) e^"" , P = b|e^" • (A.12) 

These functions satisfy the following orthonormality and completeness properties 

■(l}.AP)€'^n'iP) = 2^"" 6{V - V') 5n,n' , (A.13) 

Y. / dv^,Ap)^l^^{j/) = 4^2^(lnH2-lnbf)5(a-a')- (A.14) 

-.— .-^ J —OQ 






n— — 00 



The homogeneous BFKL equation for the octet case can be written in the form 



bP ^ 1 
7:^ + 2^ 



wV = -aHtP , F = In ^ + -Ho , (A. 15) 
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where Hq is the Hamihonian for the singlet case. Its solution can be constructed in terms of the linear 
combination of the functions 4>u,nip) 



^u^Ap) = e" 



\P?\ , . 



Taking into account (JA.ISP the function a^^A^^ should satisfy the equation 
' d , (\ \n\\ , /I |?i|\ „ ^ 



a.uj,-n(y)- 



(A.16) 



(A.17) 



(^ + 2^) 



r ( ^ + ?,;/ + -y 

-n /' 1 ■ I I" 



wa^,„(i/) = -a 
Its solution is 

auj.niy) = exp 

For tti^Ai') we have the following normalization 

/oo 
-OO 

and completeness conditions 

/oo 
dwa^A^)alA'^') ^ 2-Ka5{v-v'). 
-OQ 

Therefore the completeness relation for the eigenfunctions ipuj,7i{p) has the form 

E y ^V'.,„(p)C„b')=«(27r)2^(lnH2-ln|p'n<5(a-a'). (A.21) 



(A.18) 



(A.19) 



(A.20) 



n— — OO 



Using these expressions, we can solve the inhomogeneous equation for the Green's function gui'P^p' ^ 0) 



u;gUp,lf,0) = (2^)2 - \p\^ 5\p-p') - aH g^{p,0 ,0) 
in terms of a superposition of eigenfunctions of the homogeneous equation, i.e. 



9l 






UJ ~ UJ' 



(A.22) 



(A.23) 



where P means that the integral over uj is taken with the principal value prescription. 
In terms of this Green's function the i2-channel partial wave can be written as 



M^) = 7Z 



Z f°° d\k\^ f^'^ dak*kl [^ d\k'\^ f^'^da'k'h 



At: J, |fc|3 ^ k*+klJo \k'\ 



ko - k 



ygUk,k',0), (A.24) 



where Z is the divergent factor discussed in (jA.5[) . We can now write /2(w) in a different form using 
the explicit expression for gui{k, k' , 0): 






du' 



n— — oo 



where 



bujn{kl 



-oo 27ri(w - w') 
dvK^n{ki)a^A^) , 



/\ t>uj'n{'^) 



(A.25) 



(A.26) 
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6c^„(fc2)-/ duK^^{k2)al^^{u). (A.27) 

J — OO 

The functions K and K read 

K^n{ki)= — — ^ -^— / — — — — = i„i_i . , (A-28) 



where 






Thus, one can obtain the foUowing simple representation for the functions h and b 

bUki) - -27r(-l)"e^"i(l"l-i)c^„(fci), (A.31) 

b^nih) = -27re-^"^(l"l-i)<„(fci), (A.32) 

where 



Cu.n{ki)= I dve V" " y -1^^^^ ^. (A.33) 

2~ 

The final result for the imaginary part of the amplitude in the variable S2 reads 



^^^^^^^ = Z.Y. (-l)"e»-(l"l-i) r d^ f^y cUhXM, (A.34) 

^ n=-oo •'-°° \^ ^ 

where ai2 — ai— a2- There is an ambiguity in the integration over v at v = Q for \n\ = 1, but at that 
point Cuin{ki) does not depend on u>. 

A. 2 Spectrum quantization 

It is important to note that in the region p ^ we can not use the simplest form for the Regge 
trajectory in the Born approximation and we should write the exact expression instead: 

\n\p\^-i^E,i\p\) = -l(^^^ . (A.35) 

The reason for this is that the solution ijju:.n{p} has a good behavior only for large \p\, i.e. 

1+1" 
lim V^,„(p) ~ e^"" f ^) ' . (A.36) 

In the region of small \p\ its asymptotics is given by the saddle point contribution in the integral over 
v and is not stable for the simplified expression ii^g(|p|) at e — > 0. The position of this saddle point v 
is defined by the solution of the BFKL equation in the classical approximation 



., ,x 1 , A \n\\ (\ In, , 



(A.37) 



In this expression we have used the exact expression i?g(|p|) for the gluon energy. Let us indicate that, 
due to the symmetry oiEQ(y) under the substitution v -^ —v, there are two solutions of this equation 
related by this symmetry and the semiclassical expression for 'ipui,nip) oscillates in this region. 
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The intercept A of the corresponding singularity in the j — 1-plane of the i-channel corresponds 
to the values ly = 0,n = oi the Mobius parameters (for e < 0) 

2 

A = -amin(£;g(|p|)-41n2) = ^, ujp^ ^Ncln2, (A.38) 

where ujp is the intercept of the BFKL Pomeron. 

Let us solve the Schrodinger equation for the wave function with the modified expression for the 
Regge trajectory analytically. For this purpose we shall use the representation where the coordinate 
is 

2; = ln^. (A.39) 

The Schrodinger equation has the form 

Eoiiy, n)*,,„(x) - ('^ + ^) *.,n(a^) , (A.40) 

where e -^ —0 and Eq{v, n) is the total energy at a; — > — cx) 

i?o(^^,n)=^Q+2t^+^)+^Q-z^+^)+27. (A.41) 

At X -^ — oo the potential energy goes to zero and we can search for two solutions of this equation of 
the form 

^±„(x) = e±-- Yl ^^i") ^"''' ' (A-42) 

r=0 

where C^{v) satisfies the recurrence relation 



Eo{v,n)Ct{v) = --^Cf_M + E„{^ ±ir,n)Ct{y). (A.43) 

Therefore one can write the following expression for the coefficients C^ [v) 

^fM - {^) n g„(..„)J„(.^.,„) ■ (A^«) 

In principle, the expansion in (|A.42p with these coefficients is convergent for all values of x and 
therefore we could find at least numerically the linear combination of '5+ and ^^ for which the wave 
function decreases at a; ^ oo. 

We consider now the case of small z/, where we can use the diffusion approximation for Eq: 

i^o(.,n) = 2^(^- + Uj+2,-^"(^- + Uj_. (A.45) 

In particular for n = we have 

£:o(i/,n) = -4 1n2 + 14C(3)i'^ (A.46) 

In this case one obtains 

C±M = ( -^ \\\ ^ (-14C(3)e)-r(l±2z.) 

'■ ^ ^ \UQ{2,)e) f}^t{t±2iv) r(r + 1) r(r + 1 ± 2iz.) ' ^ ' 
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As a result, we ean express ^ with an appropriate normalization constant in terms of the Bessel 
function with imaginary argument (for e < 0) 



*j„,.,.,_Ji£|_^V ,. 



The solution, which has the good asymptotic behavior at a; ^ oo 



where 



*^,o(a;) ^ e V '(-"^asw' , (A.49) 



^M^)^K2.JJ^^^^^-^], (A.50) 



-7C(3)6 

z) - hii 
2 sin(27rw) 



K2iu{z) ^ — . , , . (A. 51) 
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B Calculation of the one and two loop contributions 

In this appendix we calculate, starting from (j93p . the one and two loop approximations. For the one 
loop approximation we compute the integral: 



OQ poo 

M= Y. i-ir / - 

, , — , J —oo ^ 



dv 



n— — oo 



where 



and 



2 ^ „2/4 

qskiqtk^ 



«'^^"/2 



P'^a 



P 



q^klqik2 



We begin with the terms n ^ and integrate over ly; 



(B.l) 

(B.2) 
(B.3) 



M, 



n^tO 



\n=l 
/ °° 9 



/2/a"/2+a-"/2 



a-/2+a-"/2 



-e{P - 1)2^ (in ((1 + \/S/V^)(l + 1/V^))) 
-0(1 - /3)27r (in ((1 + ^/V^)]{1 + ^^(3))) 



Using (B.2) and (B.3) we obtain: 



M„^o = -27r In — — ■ In /3 

* \kik2qiq3\ 2 



(B.4) 



(B.5) 



For the term n = the divergence at z^ = needs to be regularized. In order to reproduce the one 
loop result (|25p derived in dimensional regularization we need 



Regs2 



dv 

7/2 



qski 



k2qi 



2iu 



27r 



i + mH 



fi- 



In 



93 91 



fcifc2 



In 



kallfci 



Ifelkil 



(B.6) 



The sum of this contribution for n = and Af„^o is (apart from the overall factor) in agreement with 
the one loop result in (pS)) . 

Next let us consider the two-loop contribution. We need to calculate the following integral: 



R= Y. (-1)" r -?^r^/3*'^«"/' 

^_ , J — c 



n— — oo 



where according to eq. (38) and (85) 



-^un 






/-oo i^^ + n^/4: 



1 



(B.7) 



1 



fc+1 k+l + iu+\n\/2 k + l-w+\ji\/2 



(B. 



The integral over v can be calculated by residues, and we take into account that the contributions 
from the poles v = ±i|n|/2 exist only for n ^ 0. As for the other poles, they give contributions also 
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at n = 0. Thus, we obtain 

n=l \ fe=0 ^ 

n— 1 

oo oo 



, , „ x- . ^ k + l + n 

n—1 \ k—0 ^ 



oo oo _, 

.-(fc+i+M) L 



(fc + l)(fc+l + |n|) 

n— — OO 

= -7rln/3 ^ ^^ ^(a^ + a'^) {9{(3 - 1) /J^t _ 0{i _ /3)/3t) 



fc=0 

(fc + l)(fc + l + Inl) 



n 

n— 1 



-27r^(-ir(at+a-t) (-J^ + i^i) {0(0 - 1) p-^ + 9(1 ~ P)0'- 



n— 1 \ r—1 / 



(oo oo o — m °^ o- 

-7rln/3 ^ ^^ '—{a"^ + a"*) {6{I3 - 1) /J"* - e{l - P)l3^) 



"2 



ni — l n2 — 1 

+27r0(l-/3) y ^ >—a^p^ V ^^^a-^/g"? . (B.9) 

ni— 1 n2 — 1 

In obtaining two last contributions we passed to the new summation variables tii = to + ri/2 and 
772 = 771 — 77/2. These transformations give a possibility to write the total result for i? in the following 
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simple form 



R ^ 7r0(/3-l) -ln/3^^^ ^(a* + a"*)/?"* + ^ ^^ ^(a* + a^^)/?-? 



(-1) 



(-1)\ - 



+7r6'(l - /3) ln/3 ^ ^^ —{a'i + a-t)/3t + ^ ^^ ^(at + a-t)/3t 



(-1) 



(-ir 



e{l3 - 1) XI ^ ^("* + a"^)/3"^ X ^^ ^(a* + a"*)/3"* - ln/3 






)/3t j I ^ ^^(a^ + a-*)/3* + ln/3 



Vn=l 



= 7r0(/3 - 1) In ((1 + V^/v//3)(l + l/v/"/?)) (in ((1 + V^/^){1 + l/Vo^)) + In/?) 

+Tie{l - (3) In ((1 + v//3/\/^)(l + y^)) (in ((1 + ^/J3/^/^){l + V??^)) - In/?) .(B.IO) 
Using finally the above expression for M„^o we obtain 



^^ , \h + k,\'\q2\' , |fcl + fc2ng2p 
Combination of the one and two loop results leads to: 



(B.ll) 



1 






^s^(*^^ (In 



n— — oo 
1211, 12 



rfi^ ( qlkl\ 2 / qafcl \ ' „^(i^,n) 



^ Vfc2*9lV V^2gi 



fcl|'|fc2|' la, , \ki + H'\q2\' , |fci + fc2h92 



— In S2 In 



In- 



lA^i+felV e 2 IfePlgiP |fci|2|g3 



+ 0{a')]. (B.12) 



Indeed, the second term of the expansion in a coincides with the result (63), obtained by an indepen- 
dent calculation. 

In an analogous way the imaginary part in s can be written in the form 



n— — oo 



dv fqlkiX" "^ ( qzki 



/2 + f \k*2qlJ \k2q1 



iu+-. 



is. M.^, . A ,J<« Y. (-!)"«'». / V^ ( S ) " ( !^ ) " »?"""' , (B^13) 



where 



Regs I -^ 



qski 



k2qi 



2^(-i + lnM!+ln 



9391 

kik2 


- 


In 


qaki 
qik2 



(B.14) 



It corresponds to the following expansion in a 
1 



-3s M2- 



|2|„.|2 



■So I In T^ . ,0 r, In S2 In — rr~{T, — n3 — In — |^— pj^ — h 0{a 



2"^ V"|fci + fc2|V e 2 '^ |fc2P|<Zi|2 IfciPksP 

For the imaginary part of the amplitude M^^^ in the variable t'2 we obtain the similar result 



(B.15) 



1 



^%,M,^,^^t'r^^^^ j:i-irReg,, 



dv 



(m 



,,2 I "^ I h*h* 
V + ^- \ "^2 "^l 



k2ki 



9391 \ ,'c^(i^,n) 



^^o.,.,.;^ (B.16) 
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where in this case the regularized integral over v for n — Q and a — Q does not contain any 1/e 
divergence 



Regv. 



7^ 



9391 



k2ki 



2-K In 



9391 


- 


In 


9391 

fclfc2 



(B.17) 



It gives the following a-expansion of 5t' M-^^z 
1, 



-9,, M3_3 



TT 
77 t-o 



In- 



l9iPl93p _ a , , , |9i+93-9^Pl92p , |9i + 93 - 92p|92p 

2|„„|2 9^-2 11 \h.\1\h.A2 "^ l„J2U„|2 



|9i +93 -92r|92 



l^ll 



i9in93i 



o(«') 



Analogously we find the the imaginary part of the amplitude Mj,^^ in the variable s 



let A/f - " /'^(*2) 

TT 47r 



E (-l)"^eg. 



(<m 



y2^r^ yk*kl 



9391 \ ,'i^{i^,n} 

'■2 



fc2fcl 



where 



Reg, I — 



9391 



k2ki 



-2tt 



In 



9391 
kik2 


- 


In 


9391 
kik2 



(B.18) 



(B.19) 



(B.20) 



The expansion in a beyond the one loop approximation coincides with that of Sf M^^^ 



-3,M3^3 = 

t' '"'''" lln- 



|fciP|fc2| 



91 +93 -92PI92P 2 



flj^^, ^^^ l9i + g3-9^Pl92p ^^^ I9i + g3-g2p|92p 



|fe|2|A:i|2 



I91PI93P 



Oia^) 



(B.21) 
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